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BrabA.1 1339.a: anomalous diffraction and ligand 
binding guide towards the elucidation of the 
function of a 'putative /?-lactamase-like protein' 
from Brucella melitensis 

The crystal structure of a /J-lactamase-like protein from Brucella melitensis was 
initially solved by SAD phasing from an in-house data set collected on a crystal 
soaked with iodide. A high-resolution data set was collected at a synchroton at 
the Se edge wavelength, which also provided an independent source of phasing 
using a small anomalous signal from metal ions in the active site. Comparisons of 
anomalous peak heights at various wavelengths allowed the identification of the 
active-site metal ions as manganese. In the native data set a partially occupied 
GMP could be identified. When co-crystallized with AMPPNP or GMPPNP, 
clear density for the hydrolyzed analogs was observed, providing hints to the 
function of the protein. 



1 . Introduction 

The mission of the Seattle Structural Genomics Center for Infectious 
Disease (SSGCID) is to provide the scientific community with a 
library of crystal structures of proteins from infectious organisms that 
are potential drug targets. This collection of structures is hoped to 
promote structure-guided drug design against National Institutes for 
Allergy and Infectious Disease (NIAID) class A, B and C infectious 
disease organisms. The target organisms are divided between 
SSGCID and its sister organization the Center for Structural Geno- 
mics of Infectious Diseases (CSGID). Both centers provide plasmids 
(BEI repository), protein samples (SSGCID), coordinates and 
structure factors (Protein Data Bank), and raw diffraction images 
(http://www.csgid.org/csgid/pages/diffraction_images) to the commu- 
nity immediately after completion of a target. 

Brucella melitensis is a Gram-negative bacterial pathogen and one 
of the several Brucella species identified to cause brucellosis (Moreno 
& Moriyon, 2002). Brucellosis affects mostly goats, sheep and cattle, 
causing reproductive diseases such as abortions and stillbirths, hence 
the popular name 'contagious abortion'. The pathogen can be 
transmitted to humans via contact with infected animals or the con- 
sumption of unpasteurized milk or cheese (Young, 1995). In humans 
the disease was initially called 'Malta fever' or 'undulating fever' and 
is multi-systemic, with nonspecific symptoms such as fever, chills, 
dementia, fatigue, headaches, nausea, vomiting, muscular and bodily 
pain. Regions with significant impact of brucellosis are parts of the 
Mediterranean, Asia, Africa and Latin America (Corbel, 1997). 
Infections can be successfully treated with antimicrobial agents. The 
pathobiology of Brucella is unique owing to the absence of classical 
virulence factors (Moreno & Moriyon, 2002). 

Here, we report the high-resolution structure of BrabA.1 1339.a, a 
structural genomics target that entered the pipeline as a '/3-lactamase- 
like protein'. /3-Lactamases break down /3-lactam antibiotics such as 
penicillin and lead to the resistance of pathogens to antibiotics. There 
are four classes of /j-lactamase: classes A-D. While classes A, C and 
D have an active site-serine, class B metallo-/j-lactamases have a 
binuclear metal center in the active site that activates a water. With 
high-resolution data including anomalous information at various 
wavelengths, we could obtain a detailed description of the active site 
of BrabA.11339.a and an identification of the active-site metals. 
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Table 1 

Data-collection statistics. 



Values in parentheses are for the highest of 20 resolution shells. 



Ligand 


Iodide 


Native Cu Ka 


Native ALS 501 


GMP Cu Ka 


AMP Cu Ka 


AMP ALS 502 


Wavelength (A) 


1.5418 


1.5418 


0.9774 


1.5418 


1.5418 


1.3476 


Space group 


C222, 


C222, 


C222 t 


C222, 


C222! 


C222, 


Unit-cell parameters (A) 


a = 72.71, b = 74.95, 


a = 72.63, b = 75.04, 


a = 72.89, b = 75.25, 


a = 72.86, b = 75.23, 


a = 72.85, b = 75.20, 


a = 73.21, b = 75.53, 




c = 98.27 


c = 98.11 


c = 98.36 


c = 98.54 


c = 98.53 


c = 98.94 


Resolution range (A) 


50-2.00 (2.05-2.00) 


50-1.85 (1.90-1.85) 


50-1.27 (1.30-1.27) 


50-1.70 (1.74-1.70) 


50-1.70 (1.74-1.70) 


50-1.60 (1.64-1.60) 


Unique reflections 


17854 (1257) 


23011 (1520) 


71257 (5123) 


29984 (2068) 


29964 (2069) 


36362 (2574) 


Multiplicity 


12.4 (6.9) 


5.4 (1.90) 


9.5 (6.7) 


9.1 (4.4) 


9.1 (4.5) 


12.4 (6.1) 


Completeness (%) 


99.3 (94.3) 


98.9 (91.3) 


99.7 (97.2) 


99.5 (93.1) 


99.5 (93.6) 


99.6 (95.9) 


^mergef 


0.072 (0.140) 


0.048 (0.149) 


0.051 (0.489) 


0.037 (0.147) 


0.041 (0.203) 


0.069 (0.390) 


Mean IltrQ) 


27.1 (11.4) 


24.7 (5.1) 


28.3 (3.9) 


40.0 (9.2) 


35.2 (7.7) 


25.1 (5.8) 



Ligand-bound structures lead us towards an understanding of the 
function of this protein. 



2. Materials and methods 

2.1. Cloning, expression and purification 

The gene for BrabA.11339.a from B. melitensis biovar Abortus 
2308 (NCBI YP_414422.1; locus BAB1_1016; UniProt Q2YQ74) 
spanning the full-length protein from residues 2-272 was cloned into 
a pAVA0421 vector using ligation-independent cloning (LIC; 
Aslanidis & de Jong, 1990; Mehlin et al., 2006). The vector encodes an 
N-terminal histidine-affinity tag followed by the substrate sequence 
for 3C protease; the entire tag sequence is MAHHHHHHMGTLE- 
AQTQ'GPGS, where ' indicates the cleavage site. 

The protein was expressed in Escherichia coli using BL21 (DE3) 
R3 Rosetta cells and autoinduction medium (Studier, 2005) in a LEX 
Bioreactor. Starter cultures of LB broth with appropriate antibiotics 
were grown for ~18 h at 310 K. ZYP-5052 auto-induction medium 
was freshly prepared as per Studier's published protocol (Studier, 
2005). Antibiotics were added to 2 1 bottles of sterile auto-induction 
medium. The bottles were inoculated with all of the overnight culture. 
Inoculated bottles were then placed into a LEX bioreactor. Cultures 
were grown for ~24 h at 298 K; the temperature was reduced to 
288 K and the cultures were grown for a further ~60 h. To harvest, 
the medium was centrifuged at 4000g for 20 min at 277 K. The cell 
paste was flash-frozen in liquid nitrogen and stored at 193 K. 

The frozen cells were resuspended in lysis buffer (25 mM HEPES 
pH 7.0, 500 mM NaCl, 5 % glycerol, 30 mM imidazole, 0.025 % sodium 
azide, 0.5% CHAPS, 10 mM MgCl 2 , 1 mM TCEP, 250 ng ml" 1 
AEBSF and 0.05 ug ml -1 lysozyme). The resuspended cell pellet was 
disrupted on ice for 30 min with a Virtis sonicator (408912) set to 
100 W power with alternating cycles of 15 s pulse-on and 15 s pulse- 
off. The cell debris was incubated with 20 ul Benzonase nuclease 
(25 units ml -1 ) at room temperature for 45 min and clarified by cen- 
trifugation using a Sorvall SLA-1500 centrifuge at 14 000 rev min -1 
for 75 min at 277 K. The protein was purified from the clarified cell 
lysate by immobilized metal-affinity chromatography using a 5 ml 
HisTrap FF column (GE Healthcare) equilibrated with binding buffer 
(25 mM HEPES pH 7.0, 500 mM NaCl, 5% glycerol, 30 mM imida- 
zole, 0.025% sodium azide, 1 mM TCEP). The recombinant protein 
was eluted with 250 mM imidazole. The hexahistidine tag was 
removed by overnight incubation with 3C protease. The sample was 
then passed over HisTrap resin in a gravity-flow column. The tagless 
protein was collected in the flowthrough. It was further purified 
by size-exclusion chromatography (SEC) using a HiLoad 26/60 
Superdex 75 column (GE Healthcare) equilibrated with SEC buffer 



(25 mM HEPES pH 7.0, 500 mM NaCl, 2 mM DTT, 0.025% sodium 
azide and 5% glycerol). The protein eluted as a single peak. The peak 
fractions were pooled and concentrated to 19.6 mg ml -1 using an 
Amicon Ultra centrifugal filter with a 10 kDa cutoff (Millipore). The 
concentrated sample was flash-frozen in liquid nitrogen and stored at 
193 K. 



2.2. Crystallization and X-ray data collection 

For the target BrabA.11339.a at a protein concentration of 
19.6 mg ml -1 , two sparse-matrix screens were set up following the 
strategy of Newman et al. (2005): JCSG+ (Emerald BioSystems) and 
PACT (Molecular Dimensions). 0.4 pi protein solution was mixed 
with 0.4 pi well solution and equilibrated against 100 pi reservoir 
using 96-well Compact Jr crystallization plates (Emerald Bio- 
Systems). Crystals suitable for diffraction studies were obtained 
under two very similar conditions from the PACT screen: G9 
(200 mM sodium/potassium tartrate, 100 mM Bis-Tris propane pH 
7.5, 20% PEG 3350) and H9 (200 mM sodium/potassium tartrate, 
100 mM Bis-Tris propane pH 8.5, 20% PEG 3350). Crystals were 
cryoprotected in a buffer consisting of well solution to which 20% 
ethylene glycol was added. 

A native data set was collected in-house at Cu Ka wavelength 
using a Rigaku SuperBright FR-E + rotating-anode X-ray generator 
equipped with Osmic VariMax HF optics and a Saturn 944+ CCD 
detector (Table 1, native Cu Ka). The diffraction data were reduced 
with XDSIXSCALE (Kabsch, 2010) and scaled to 1.85 A resolution 
in space group C222 t . 

In order to obtain experimental phase information, a crystal from 
the same well was incubated for 5 min in an iodide-containing buffer 
consisting of 200 mM sodium/potassium tartrate, 100 mM HEPES pH 
7.0, 25% PEG 3350, 250 mM KI. A data set with high multiplicity was 
collected in-house using the same set up as described above. 720 
images were collected with a width per frame of 0.5°. The diffraction 
data were reduced with XDSIXSCALE (Kabsch, 2010) and scaled 
to 2.0 A resolution. Data analysis in XSCALE (Anomalous Corre- 
lation and SigAno) indicated a strong anomalous signal out to 2.0 A 
resolution. 

A high-resolution data set was collected on beamline 5.0.1 of the 
Berkeley Center for Structural Biology in the context of the Colla- 
borative Crystallography program. The wavelength of this beamline 
is optimized for an Se-SAD experiment: 0.9774 A/12.4 keV (Morton 
et al., 2007). The beamline is equipped with a 3 x 3 tiled ADSC Q315r 
detector. The diffraction data were reduced with XDSIXSCALE 
(Kabsch, 2010) and scaled to 1.27 A resolution. A weak anomalous 
signal can be detected in the XSCALE data analysis. 
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Table 2 

Refinement and model statistics. 



The models for the iodide and native Cu Ka data sets have only been partially refined. Preliminary refinement statistics are shown in italics. Values in parentheses are for the highest of 20 
resolution shells. 



Ligand 


Iodide 


Native Cu Ka 


Native ALS 501 


GMP Cu Ka 


AMP Cu Ka 


AMP ALS 502 


Resolution range (A) 


50-2.00 (2.05-2.00) 


50-1.85 (1.90-1.85) 


50-1.27 (1.30-1.27) 


50-1.70 (1.74-1.70) 


50-1.70 (1.74-1.70) 


50-1.60 (1.64-1.60) 


^cryslt 


0.242 (0.226) 


0.181 (0.213) 


0.107 (0.152) 


0.140 (0.163) 


0.164 (0.174) 


0.144 (0.223) 


Afreet 


0.295 (0.281) 


0.228 (0.266) 


0.134 (0.203) 


0.168 (0.206) 


0.196 (0.211) 


0.168 (0.260) 


R.m.s.d. bonds (A) 


0.017 


0.023 


0.018 


0.013 


0.015 


0.015 


R.m.s.d. angles (°) 


1.65 


1.94 


1.87 


1.52 


1.64 


1.65 


Protein atoms 






2116 


2116 


2110 


2118 


Nonprotein atoms 






380 


381 


344 


344 


Mean B factor (A 2 ) 






9.85 


8.13 


10.4 


10.6 


Residues in favored region 






243 [98%] 


251 [97%] 


248 [96%] 


253 [98%] 


Residues in allowed region 






4 [1.4%] 


6 [2.3%] 


9 [3.5%] 


5 [1.9%] 


Residues in disallowed region 






1 [0.4%] 


1 [0.4%] 


1 [0.4%] 


1 [0.4%] 


MulPrubity score [percentile] 






1.92 [92nd] 


1.27 [98th] 


1.38 [96th] 


1.08 [99th] 


PDB code 


n/a 


n/a 


3md7 


3py6 


3py5 


3qh8 



t Rc, 



2^1hki [l^obsl ~ l^caicl |/X)a*z l^obsl- The free R factor was calculated with an equivalent equation using 5% of the reflections that were omitted from the refinement. 



In order to validate ligand binding, the protein was incubated with 
1 mM AMPPNP or GMPPNP and 2 mM MgCl 2 . A small optimization 
screen was set up: 15-20% PEG 3350, 200 mM sodium/potassium 
tartrate, 100 mM Bis-Tris pH 7.0-8.0. After crystals had grown to full 
size within a few days, they were incubated for 1 min in a cryopro- 
tection buffer consisting of 1 mM AMPPNP or GMPPNP 2 mM 
MgCl 2 , 15% ethylene glycol, 21% PEG 3350, 165 mM sodium/ 
potassium tartrate, 83 mM Bis-Tris pH 7.5. For each of the co-crystals, 
a data set with high multiplicity was collected at the home source as 
described for the iodide soak above. Each data set was scaled to 1.7 A 
resolution. For another AMPPNP cocrystal from the same optimi- 
zation screen, a data set with high multiplicity was collected at a 
wavelength of 1.3476 A on beamline 5.0.2 of the Berkeley Center for 
Structural Biology in the context of the Collaborative Crystal- 
lography program. The wavelength was selected to be at lower energy 
than the Zn edge (1.284 A/9.659 keV) but at significantly higher 
energy than the Mn edge (1.896 A/6.359 keV) The ALS beamline is 
equipped with a 3 x 3 tiled ADSC Q315r detector. The data were 
scaled to 1.6 A resolution. 

2.3. Structure solution and refinement 

The packing density of the crystals (Matthews, 1968), 
Vjvi = 2.2 A 3 Da -1 (44% solvent content), suggested one molecule of 
BrabA.11339.a per asymmetric unit. A search of the Protein Data 
Bank with BLAST (Altschul et al, 1997) did not yield any structural 
homologs beyond 23% sequence identity. Attempts to solve the 
structure of BrabA.11339.a by molecular replacement with Phaser 
(McCoy et al, 2007) using the initial native data set (native Cu Ka, 
1.85 A resolution) and PDB entry 3glp (23% identity; K. Podze- 
linska, S. He, M. Wathier, A. Yakunin, M. Proudfoot, B. Ffove-Jensen, 
D. Zechel & Z. Jia, unpublished work) as search models were not 
successful. Experimental phases were obtained independently from 
both the iodide soak and the high-resolution synchrotron data set. 

2.3.1. Iodide data set. Data analysis with XSCALE indicated a 
strong anomalous signal to 2.0 A resolution, phenix.hyss (Grosse- 
Kunstleve & Adams, 2003) was used to locate 11 anomalous 
scatterers in the data set from the iodide soak. The anomalous sub- 
structure was refined with Phaser (McCoy et al, 2007), which then 
extended the anomalous substructure to 16 sites and yielded a figure 
of merit (FOM) of 0.49 to 2.0 A resolution. The experimental phases 
obtained from Phaser were improved with the CCP4 program Parrot 
(Cowtan, 2010). An initial model (256 residues, i\ work = 0.307) was 
built with ARPIwARP (Langer et al, 2008), followed by iterative 



rounds of manual real-space fitting in Coot (Emsley et al, 2010) and 
reciprocal-space refinement with REFMAC5 (Murshudov et al, 
2011). 

2.3.2. Native synchrotron data set. In a separate approach, the 
native high-resolution synchrotron data set was used to phase the 
BrabA.11339.a structure. The XSCALE analysis had indicated a 
moderate anomalous signal to ~2.0 A resolution, phenix.hyss 
(Grosse-Kunstleve & Adams, 2003) was able to locate two strong 
anomalous scatterers when the diffraction data were limited to 3 A 
resolution. Phaser (McCoy et al, 2007) refined and extended the 
anomalous substructure to three atoms and yielded an FOM of 0.360 
using data to 1.8 A resolution. After density improvement with 
Parrot (Cowtan, 2010), ARPIwARP (Langer et al, 2008) could build a 
good initial model with 267 residues and a R„ OTk of 0.200 with data to 
1.7 A resolution. 

The model for the high-resolution structure of BrabA.11339.a 
contains one protein molecule spanning residues Ser3-Val272, 
two manganese ions, a potassium ion, a sodium ion, a guanosine 
5-monophosphate (GMP), a tartrate molecule and 348 water mole- 
cules. The final structure was validated with Coot (Cowtan, 2010) and 
MolProbity (Chen et al, 2010). The structure was deposited in the 
PDB with the identifier 3md7 (see Table 2). 

The three structures cocrystallized with AMPPNP or GMPPNP 
were solved by direct refinement of the high-resolution structure 
from which all nonprotein atoms had been removed. Data and 
refinement statistics are summarized in Tables 1 and 2. 



3. Results and discussion 

The target BrabA.11339.a entered the SSGCID pipeline as a 
'/3-lactamase-like' protein from B. melitensis. A search of sequence 
databases for homologies yielded a variety of annotations, primarily 
other '/Mactamase-like' proteins, metal-dependent hydrolases, un- 
characterized proteins, PhnP (a phosphodiesterase of the carbon- 
phosphorus lyase pathway for phosphonate degradation) proteins 
and tRNase Z proteins. However, the sequence is highly conserved 
between various species, with many proteins sharing more than 50% 
sequence identity. A search of the Protein Data Base yielded a 
striking lack of structures with high sequence homology. The phos- 
phodiesterase of the carbon-phosphorus lyase pathway (PhnP; PDB 
entry 3glp) from E. coli is the closest homologue; PhnP shares 23% 
sequence identity with BrabA.11339.a, with a large number of inserts 
and deletions. 
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BrabA.11339.a readily crystallized in a well diffracting C-centered 
orthorhombic crystal form. The structure could not be solved by 
molecular replacement. SSGCID's first approach to obtaining 
experimental phases is quick soaks with iodide (Abendroth et al, 
20f0, 201f ), following the ideas of Nagem and Dauter (Nagem et al, 
200f). The principles behind this approach are (i) the strong anom- 
alous signal of iodide at Cu Ka X-ray wavelength (f" = 6.8; Bricogne 
et al., 2003), (ii) the compatibility of iodide with weak binding sites 
of different natures, both hydrophobic and basic, and (iii) the high 
concentrations of iodide that a crystal lattice can tolerate. In the 
SSGCID context, this experimental phasing strategy has been 
extremely successful (Abendroth et al., 2010, 2011). For 




Figure 1 

Structure of BrabA.11339.a dimer. In the ribbon representation of the crystallo- 
graphic dimer the /^-strands for each protomer are colored slightly differently. The 
crystallographic dyad runs vertically. Various ions are shown as spheres: Mn 2+ in 
magenta, K + in blue and Na + in green. Five strong iodide sites are shown as orange 
spheres for one protomer. GMP and tartrate are shown as stick models. 



BrabA.11339.a the phasing information was not only obtained from 
iodide ions: two metal ions in the active site and a potassium ion 
contributed in addition to the iodide ions. In fact, the high-resolution 
native data set was collected on a beamline which was optimized for 
data collection at the Se edge (0.9774 A). The minor anomalous 
signal from these three anomalous scatterers was strong enough to 
phase the structure independently. The structure solution of 
BrabA.11339.a underlines the power of anomalous phasing by both 
intrinsic anomalous scatterers such as metals and externally added 
scatterers such as iodide. 

BrabA.11339.a crystallizes with one monomer per asymmetric unit. 
A dimer is formed by crystallographic symmetry (Fig. 1). The dimer 
buries 1444 A 2 accessible surface per protomer in the interface 
as calculated by PISA (Krissinel & Henrick, 1997). BrabA.11339.a 
assumes a metallo-yS-lactamase-like fold with two central /J-sheets 
(Fig. 2). A search for structural homology using SSM (Krissinel & 
Henrick, 2004) reveals as the closest structural homologs the PhnP 
protein (PDB entry 3glp, 1.7 A r.m.s.d., 25% sequence identity), a 
putative ribonuclease (PDB entry lzkp, 1.9 A r.m.s.d., 21% sequence 
identity; Midwest Center for Structural Genomics, unpublished work) 
and ribonuclease Z (PDB entry ly44, 2.3 A r.m.s.d., 21% sequence 
identity; de la Sierra-Gallay et al., 2005) (see Fig. 2). The structural 
homology is focused mostly around the central /S-sheets and the 
helices a4 and a5, while there is more structural variety between 
,6-strands ,65 and ,66. 

As a very typical feature of /3-lactamase-like proteins and the 
phosphodiesterase protein family, a two-metal center is present in the 
active site of BrabA.11339.a (Fig. 3). The two metals are 3.40 A apart, 
bridged by water 330 and liganded by several histidine residues 
(His86, His88, Hisl70 and His241) with bond distances between 2.22 
and 2.24 A, aspartate residues (Asp90 and Aspl88) with bond 
distances between 2.28 and 2.31 A and the phosphate group of GMP. 





3g1p 




' 1zkp 



1y44 



BrabA.11339.a 



Figure 2 

Fold comparison. Comparison of the fold of BrabA.11339.a with its structural homologues. 3glp is the PhnP protein from E. coli, lzkp is a putative ribonuclease from 
B. anthracis and ly44 is ribonuclease Z from B. subtilis. The structural homology is focused around the two ^-sheets and helices a4 and oe5. Color scheme: Mn, magenta; Zn, 
cyan; Na, green; K, blue 
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Table 3 

Bond distances of the dinuclear Mn 2+ center for the high-resolution data set. 



Atom 


Distance from Mn300 (A) 


Atom 


Distance from Mn301 (A) 


Mn301 


3.40 


Mn300 


3.40 


Wat330 


2.1(1 


Wat330 


2.07 


GMP OIP 


2.38 


GMP 03P 


2.20 


Asp90 OD2 


2.28 






Aspl88 OD2 


2.28 


Aspl88 OD2 


2.31 


Aspl88 OD1 


3.00 






His24f NE2 


2.23 


His86 NE2 


2.24 






His88 ND1 


2.22 






Hisl70 NE2 


2.22 



The bond distances around the metal center are summarized in 
Table 3. They fall within the distances expected for the coordination 
of manganese ions (Zheng et ah, 2008). 

Podzelinska et ah (2009) have performed a detailed functional 
analysis of the metal-dependency of PhnP; they report a strong 
preference for Mn 2+ and Ni 2+ (Podzelinska et al. , 2009) and model the 
metals as manganese ions in their structure. In the lzkp structure of 
a putative ribonuclease (Midwest Center for Structural Genomics, 
unpublsihed work) the metals are modeled as Zn 2+ ions, likely 
because of the presence of zinc chloride in the crystallization buffer. 
In the structure of RNAse Z (PDB entry ly44) the metals are 
modeled as Zn 2+ ions because of functional studies, even though this 
is not further discussed in the structure paper and the B factors of the 
metal ions are higher than for the protein environment (de la Sierra- 
Gallay et al., 2005). 

The iodide data set already indicated significant anomalous density 
for the two active-site metals and the potassium site. In order to 
analyze the metal center in more detail, diffraction data sets with high 
multiplicity were collected at long X-ray wavelengths from the ligand 
co-crystals: in-house at A = 1.5418 A and at the synchrotron at 
A. = 1.3476 A. Anomalous Fourier peak heights for the metal sites are 
summarized in Table 4(a) and are compared with the theoretical /" 
coefficients for various atoms in Table 4(b). Potassium was part of the 
crystallization, soaking and cryobuffer at high concentrations that 
were consistent between the data sets. The occupancy of the potas- 
sium ion can therefore assumed to be constant and can be considered 
as a reference. Interestingly, the potassium-ion site is not conserved 
between the structural homologs. For both long-wavelength data sets 




Figure 3 

Dinuclear manganese center. The active site of BrabA.11339.a consists of two 
manganese ions (magenta), which are bridged by a water molecule (red). The 
phosphate group of GMP is in close proximity to this activated water molecule. The 
cr A -weighted 2F a — F c electron density for the high-resolution data set is contoured 
at la (blue) and the corresponding F a — F c electron density is contoured at ±3cr 
(green/red). The anomalous Fourier electron density is contoured at 10o~ (orange). 
The strong anomalous Fourier density at all three wavelengths indentifies the 
active-site metal as manganese. For bond distances refer to Table 3; for anomalous 
peak heights refer to Table 4. 



Table 4 

Anomalous Fourier peak heights and Fourier coefficients. 

(a) Anomalous Fourier maps were calculated over the full resolution range for all three 
data sets. Peaks were found with the CCP4 program PEAKMAX. 



Atom 


0.9774 A 


1.3476 A 


Cu Ka AMP 


Cu Ka GMP 


Mn300 


33.8 


31.8 


36.6 


45.8 


Mn301 


42.4 


42.1 


40.5 


47.2 


K302 


8.8 


12.6 


11.7 


14.3 


Cysl5 SG 


6.4 


5.7 


7.2 


9.0 


Cys31 SG 


6.5 


4.9 


8.3 


10.0 


Met71 SD 


5.7 


8.3 


9.2 


9.9 


Metl09 SD 


6.7 


6.6 


7.3 


7.9 


Cysl28 SG 


5.3 


5.1 


8.0 


10.3 


Metl43 SD 


5.0 


6.0 


5.7 


5.4 


Cysl86 SG 


6.5 


5.5 


7.5 


9.3 


Met242 SD 


6.1 


7.3 


8.3 


10.4 


Met252 SD 






5.9 


5.6 


(b) Anomalous Fourier coefficients /" (e ) for selected elements at Cu Ka wavelength 


and 12.4 keV are listed for comparison (Bricoj 


me et at, 2003). Manganese is the only 


metal in question that has a stronger signal than potassium at both 


wavelengths. 




0.9774 A 




1.3476 A 


1.5418 A 


I 


3.18 




5.49 


6.83 


S 


0.23 




0.43 


0.56 


Mil 


1 .30 




2.26 


2.80 


Zn 


2.48 




0.55 


0.71 


K 


0.46 




0.84 


1.08 


Ni 


1 .95 




3.30 


0.51 



and the Se-edge data set strong anomalous Fourier peaks can be 
observed for the two active-site metals. Each of the peaks is about 
four times as strong as the signal for potassium. Manganese is the 
only metal in question that has a significant anomalous signal at all 
three wavelengths. The anomalous signal for Zn 2+ is strong at the 
Se edge, but is very weak at either of the longer wavelengths. The 
anomalous signal for Ni 2+ is strong at the Se edge and at 1.3476 A; 
however, it is low at 1.5418 A. The bond distances between the metals 
are in the range expected for Mn 2+ ions (Zheng et al., 2008). There- 
fore, the dinuclear center of BrabA.11339.a can be assigned as two 
manganese ions with a high degree of confidence. 

In the native structure, a partially occupied GMP molecule can be 
observed in the immediate proximity of the binuclear Mn 2+ center. At 
no stage during the purification was GMP added to the protein; thus, 
it must have been copurified from the expression host. In order to 
confirm the binding of nucleotides, crystals were grown in the 
presence of GMPPNP and AMPPNP. In both cases, however, very 
solid electron density could be observed for GMP and AMP, 
respectively (see Fig. 4). The a-phosphate group is very well defined, 
while there is no density for the /i-phosphate or /-phosphate groups. 
The /-phosphate groups of GMPPNP or AMPPNP, which had been 
added during crystal growth and had been supplemented for a short 
soak in the cryobuffer, are nonhydrolyzable; the /i-phosphate groups, 
however, can be hydrolyzed. It can be assumed that BrabA.11339.a is 
able to hydrolyze AMPPNP and GMPPNP between the a-phosphate 
and the /S-phosphate to AMP or GMP, respectively. 

The GMP or AMP molecules are bound in the dimer interface and 
make multiple interactions. The phosphate group is involved in 
hydrogen bonds to both manganese ions, the bridging water molecule 
330 (called water 1 in the ligand-bound structures), histidine residues 
His88 and His241, aspartate residue Asp90 and further water mole- 
cules (315, 454 and 511 in the high-resolution structure). The ribose 
group makes two interactions with the carboxyl group of Asp96 (02* 
and 03*) and two interactions with the guanidine group of Arg99 
(02*), each from the other protomer. The guanine group stacks with 
Tyrl27 from the other protomer. The guanine moiety engages in 
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Figure 4 

AMP-bound and GMP-bound structures. BrabA.11339.a was cocrystallized with AMPPNP and GMPPNP, variants of AMP or GMP that are nonhydrolyzable between the 
/J-phosphate and y-phosphate. Left panels, OMfT densities for the AMPPNP (top) and GMPPNP (bottom) cocrystals (see Fig. 3 for colors). The tr A -weighted 2F 0 — F c 
electron density for the high-resolution data set is contoured at ley (blue) and the corresponding F 0 — F c electron density is contoured at ±3cr (green/red). The right panels 
show the refined densities for the models refined with AMP (top) and GMP (bottom), respectively. No density is visible beyond the a-phosphate. It is likely that 
BrabA.11339.a has hydrolyzed the /3-phosphate. The nucleotides are bound in the dimer interface and engage in multiple interactions with both protomers. 



hydrogen bonds to the amide N atom of Ala89 (N3) and to water 
molecules. Neither of these interactions should be able to provide 
specificity for either AMP or GMP. 

In the PhnP structure (PDB entry 3glp), a malonate is modelled in 
the space that is occupied by the phosphate and the ribose. No ligand 
was modelled in the putative ribonuclease (PDB entry lzkp). In the 
structure of ribonuclease Z from B. subtilis (PDB entry ly44, r.m.s.d. 
2.3-2.4 A, 20% sequence identity) a phosphate ion is modelled in an 
identical location to the phosphate group of GMP in BrabA.11339.a, 
even though the B factors for the phosphate ion are significantly 
higher than those for the environment (de la Sierra-Gallay et ai, 
2005). None of these three homologous structures are able to 
accommodate a GMP or an AMP molecule in the conformation 
observed in BrabA.11339.a without conformational changes. 

The discovery of a purine monophosphate in the vicinity of the 
dinuclear manganese center allows speculation regarding the func- 
tion of BrabA.f 1339. a. For various enzymes with a dinuclear center it 
is known that the bridging water is activated by the two metal ions. 
The water reacts as a hydroxide, starting a hydrolytic reaction. The 
phosphate group of the GMP or AMP is located in close proximity to 
this deprotonated water molecule 303 (2.2 A). There is enough space 
in the active site to fit two more phosphate groups. It therefore seems 
likely that BrabA.11339.a is a GTP/ATP or GDP/ADP hydrolase. 

The genetic context of BrabA.11339.a (BAB1_1016) hints at a 
function in nucleotide processing. Its genetic neighbors are 
BAB1_1012, DNA polymerase II subunit A; BAB1_1013, a hypo- 
thetical protein; BAB1_1014, metG, a methionyl-tRNA synthetase; 
BAB_1015, a TatD-related deoxyribonuclease; BAB1_1016, 
BrabA.11339.a; BAB1-1017, a cyclic ,2-glucan ABC transporter; 
BAB1J018, a hypothetical protein; and BAB1J019, an RNA- 
binding S4:pseudouridine synthase. 

The structure of BrabA.11339.a and its binuclear manganese 
center were determined with high accuracy using high-resolution data 



and anomalous diffraction. The initially serendipitous discovery of 
GMP in the active site and the subsequent high-resolution co-crystal 
structures with in situ hydrolyzed AMP or GMP shed light on the 
function of this protein, which entered the SSGCID pipeline as a 
'/Mactamase-like protein'. Further analysis will be needed to validate 
the biological function of this protein and its suitability for structure- 
guided drug design against the pathogen B. melitensis. 
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